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We study the stabihty of a two-component Bose-Einstein condensate (BEC) in the parameter 
regime in which its classical counterpart has regular motion. The stability is characterized by the 
fidelity for both the same and different initial states. We study as initial states the Fock states with 
definite numbers of atoms in each component of the BEC. It is found that for some initial times the 
two Fock states with all the atoms in the same component of the BEC are stabler than Fock states 
with atoms distributed in the two components. An experimental scheme is discussed, in which the 
fidelity can be measured in a direct way. 
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I. INTRODUCTION 

In many research fields, such as Bose-Einstein conden- 
sation (BEC) and quantum information processing, sta- 
ble and coherent manipulation of quantum states is of 
crucial importance [1(1. In fact, the instability issue of 
BEC in dilute gases Q has been constantly addressed 
for its crucial role in the control, manipulation, and even 
future application of this newly formed matter, including 
dynamical instability [3j] , Landau or superfluid instability 
[j| , modulation instability [1, Q , and quantum fluctua- 
tion instability . It is found that instability may break 
the coherence among the atoms and lead to collapse of 
BEC 0. 

Recently, the stability of the quantum motion of BEC 
systems under small perturbation, which is measured by 
the so-called quantum Loschmidt echo or fidelity, has 
been studied [a, Here, the fidelity is defined as 

the overlap of two states obtained by evolving the same 
initial state under two slightly different Hamiltonians 
[ni,[ll,[ll- Explicitly, it is M{t) = |m(t)p, where m(t) 
is the fidelity amplitude for an initial state \^o), defined 
as 

m(t) = {^o\exp{iHt/h)ex.p{-iHQt/h)\^Q). (1) 

Here Hq and H = Hq + eV are the unperturbed and per- 
turbed Hamiltonians, respectively, with e a small quan- 
tity and V a generic perturbing potential. 

Fidelity decay has been well studied in quantum sys- 
tems whose classical counterparts have chaotic motion 
[11 [11 [11, [13, [S [11 [13, [El, [11 m. Related to the 
perturbation strength, previous investigations show the 
existence of at least three regimes of fidelity decay, (i) 
In the perturbative regime with sufficiently weak pertur- 
bation, in which the typical transition matrix element is 
smaller than the mean level spacing, the fidelity has a 
Gaussian decay [H [H, [l3|- (h) Above the perturba- 



tive regime, the fidelity has an exponential decay with a 
rate proportional to e^, usually called the Fermi-golden- 
rule (FGR) decay of fidelity [ij, [H, [H, [13, [S [ll ■ (iii) 
Above the FGR regime is the Lyapunov regime in which 
M{t) usually has an approximate exponential decay with 
a perturbation-independent rate; the decay rate of the 
fidelity is given by the Lyapunov exponent of the under- 
lying classical dynamics, when the classical counterpart 
of the quantum system has a homogeneous phase space 
J_4, 15, 16, 19, 20, 21]. 

For quantum systems whose classical counterparts 
have regular motion, many investigations in the decaying 
behavior of fidelity have also been carried out (see, e.g., 
Ref. I, [S [H [13, [13, [H, [H, [13, [H, [H, [13, [Si), how- 
ever, the situation is still not as clear as in the case of 
quantum chaotic systems. The point is that fidelity decay 
in quantum regular systems exhibits notable initial-state 
dependence. (In quantum chaotic systems, the main fea- 
ture of fidelity decay is initial-state independent beyond 
a short initial time.) The most thoroughly studied ini- 
tial states are narrow Gaussian wave packets (coherent 
states) , for which the semiclassical theory and numerical 
simulations show that the fidelity has, roughly speaking, 
a Gaussian decay followed by a long-time power law de- 
cay [H, 113, HH ■ Some other types of initial states of prac- 
tical interest have also been studied numerically, e.g., the 
fidelity of an initial maximally entangled (N-GHZ) state 
is shown to have an interesting oscillating behavior [^. 

The above considerations motivate our interest in the 
stability of a two-component BEC (3^ , which is exposed 
to a pulsed laser field coupling two internal states of the 
atoms in the BEC. This BEC system possesses a clas- 
sical counterpart, which has chaotic or regular motion 
depending on both the strength of the coiipling field and 
that of the interaction among the atoms In Ref. [§], 
the fidelity of initial coherent states in this BEC system 
has been studied and found in agreement with previous 



2 



analytical predictions in both cases with chaotic and reg- 
ular motion in the classical limit. We remark that, more 
recently, the fidelity approach has also been employed in 
the study of the stability of another BEC system, which 
starts from the ground state [lol |. 

In this paper, we study the stability of the quantum 
motion of the same BEC system as in Ref. Q- How- 
ever, here we are interested in the stability of initial Pock 
states, which have definite numbers of atoms occupying 
each of the two components of the BEC, and in the pa- 
rameter regime in which the classical counterpart of the 
BEC system has regular motion. In particular, we are 
interested in whether some Fock states are more stable 
than others. Knowledge about the stability properties of 
this type of initial states may be useful for potential ap- 
plication of BEC in fields such as quantum information 
m. As in Ref. |], we still consider the stability issue 
caused by small perturbation due to imperfect control of 
the coupling field and use fidelity to characterize the sta- 
bility, with the difference between H and Hq in Eq. ([T]) 
given by a small change in the strength of the coupling 
field. 

Presently, there is no analytical prediction for fidelity 
decay of initial Fock states, therefore, our investigation 
is mainly based on numerical simulations. Interestingly, 
our numerical results show that, for some initial times, 
initial Fock states with all the atoms in the same compo- 
nent of the BEC are more stable than other Pock states. 
To see whether or not this property is due to the specific 
measure given in Eq. ([T]), we have also studied the be- 
havior of a more general form of the fidelity, which is the 
overlap of the time evolution of two different initial Pock 
states under two slightly different Hamiltonians. Our 
numerical simulations for this more general fidelity give 
consistent results. 

The paper is organized as follows: In the second sec- 
tion, we discuss briefly the two-component BEC model. 
Section III is devoted to a study of fidelity for the same 
initial Fock states. In Sec. IV, we introduce and study 
numerically the more general fidelity mentioned above 
for two different initial states. Conclusions are given in 
Sec. V. An experimental scheme for measuring fidelity 
of the quantum motion of a two-component BEC is dis- 
cussed in the Appendix. 



II. PHYSICAL MODEL 

We consider the same BEC system as in Ref. [§| , specif- 
ically, cooled ^^Rb atoms with two different hyperfine 
states F — l,mF — —1 and F = 2,niF = +1. The 
total number of the atoms in the BEC is N. A near res- 
onant, pulsed radiation laser field is used to couple the 
two internal states. Within the standard rotating-wave 
approximation, the Hamiltonian describing the transition 




FIG. 1: Decay of fidelity for K ^ 1, = 0.2, L = 100, and 
a = 0.1. The initial states are Fock states \l) of I = —100 (thin 
solid curve), -75 (dashed curve), (dotted curve), 75 (dashed- 
dotted curve), and 100 (thick solid curve), respectively. The 
initial decay of the fidelity of the two states \L) and | — L) is 
much slower than that of the other three Fock states. 

between the two internal states reads 

= ^(aloi - 0202) + -[a[ai - 0^02) 

+ :|(ata2 + 4Si)^5(<-nT), (2) 

n 

where K is the coupling strength proportional to the laser 
field. Here we suppose that the laser field used to couple 
the two states is turned on only at certain times with 
a period T. The operators di, d|, 0,2, and dj are boson 
annihilation and creation operators for the two compo- 
nents, respectively. The parameters are K = hriR, g — 

^2^77(2012-011-022), M = -(^+(4A^7r?i^/m)77(oii-022)- 
Here, flji is the Rabi frequency; o^j is the s-wave scat- 
tering amplitude; S is the detuning of lasers from res- 
onance, very small and negligible in our case; m is the 
mass of atom; 77 is a constant of order one independent of 
the hyperfine index, relatingto an integral of equilibrium 
condensate wave function [3J| . 

The above Hamiltonian can be written in terms of the 
SU(2) generators [H, 

- 0^0,2 -I- 0-201 0^0,2 — djdi a\ai — 0302 

" " 2 ' ^ " 2i ' " " 2 ' 

(3) 

which gives 

H = ^iL,+gLl + KSTit)L,. (4) 

The Ploquet operator describing the quantum evolution 
m one period is [i, ^ 

U = exp[-z(/^L, + gLl)T] cxp(-^i^L,), (5) 
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FIG. 2: Decay of fidelity for K = l,gc ^ 0.2, L = 100, and 
cr = 0.1. Tlie initial states are Fock states |/) of I = 100, 99, 98 
and 97, respectively. 



where the Planck constant is set unit unless otherwise 
stated. Since the overall scaling of the Hamiltonian does 
not influence dynamical properties of the system, we will 
set n unit. 

The Hilbert space for the system is spanned by the 
eigenstates of L^, denoted by |/) with I — —L,—L + 
1,...,L, where L — N/2. These states \l) are the 
Fock states. Using Ni and A''2 to denote the numbers 
of the atoms in the two components, respectively, with 
N = N1+N2, we have {Ni- N2)/2. Hence, \l) is the 
state with iVi = (L + l) atoms in the first component and 
with iV2 — {L ~ I) atoms in the second component. In 
particular, the two states with all the atoms in one of the 
two components are | — L) and \L). The SU(2) represen- 
tation of the system discussed above is quite convenient 
for the study of properties of the Fock states \l). 

Note that for some specific choice of the parameters 
the system degenerates to the quantum kicked top model 
[3^ . As in the kicked top model, an effective Planck 
constant can be introduced, TiaS = 1/^, which will be 
written as h in what follows for brevity. The system has 
a classical counterpart in the limit TV — > 00. 



III. FIDELITY DECAY FOR INITIAL FOCK 
STATES 

As mentioned in the Introduction, we consider a small 
perturbation due to imperfect control of the coupling 
field. Specifically, for an unperturbed Hamiltonian Hq 
with the form given on the right-hand side of Eq. ^ , the 
perturbed Hamiltonian {Hq + eV) is given by the change 
K ^ K -\- e. Denoting the one-period evolution opera- 
tors corresponding to the unperturbed and the perturbed 
Hamiltonians by U and [/e, respectively, the fidelity am- 




FIG. 3: Values of M{t) for t = 1000, as a function of K. 
The five curves correspond to five initial Fock states: the thin 
solid curve for \ — L), dashed curve for | — L/2), dotted curve 
for |0), dashed-dotted curve for \L/2), and thick sold curve 
for \L). Parameters are L = 100, gc = 0.2, and a — 0.01 for 
the upper panel and a = 0.04 for the lower panel. 

plitude m(t) is now written as 

m(t = nT) = ($o| {Uiy o |$o), (6) 

with |<I>o) indicating the initial state. Fast decay of the 
fidelity means rapid lose of information during the quan- 
tum evolution in the presence of the perturbation. For 
small e, it is usually convenient to use a = e/h as a mea- 
sure for the strength of quantum perturbation. 

In this paper, we consider only the parameter regime 
in which the corresponding classical system has regular 
motion (see Ref. [9] for details of the regime). We have 
carried out numerical investigations in the fidelity decay 
of initial Fock states \l) in this parameter regime. It is 
found that the two Fock states with |/| = L, i.e., with all 
the atoms occupying the same component of the BEC, 
behave differently from other Fock states. 

Some examples of our simulations are shown in Fig. [U 
with parameters K = 1, gc = gL = 0.2, L = 100, and 
a = 0.1. For some initial times, the fidelity M{t) of the 
two initial Fock states with |/| = L has a decay which is 
approximately a Gaussian decay and is much slower than 
the fidelity decay of the other three Fock states. For long 
times, the fidelity of \L) and \ —L) are smaller than that of 
the other three Fock states, the latter of which oscillates 
and decays slowly on average. Similar results have also 
been found for other Fock states \l) with \l \ not close to L. 
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herent states is known to have initial Gaussian decay [I^ . 
Hence, it is natural to check the relationship between the 
Fock states and the coherent states generated by gener- 
ators of the group SU(2). A SU(2) coherent state \a) 
centered at a point in the sphere with polar angle 6 and 
azimuthal angle is given by l37i, jiSj 



— L), with a = 



Tr-0 



(7) 



The state | — L) is a coherent state with a = 0. To see 
the relation of other Fock states to coherent states, one 
can expand \a) in the Fock states 



*\l+L 



(1 + zz*) 



(2L)! 



[L + iy.iL-iy. 



1/2 



10, (8) 



FIG. 4: Variation of |(q|/>P with 61 for L = 100, where |a) 
is the coherent state in Eq. For the three curves on the 
left, I = 100 (solid curve), 99 (dashed curve) and 98 (dotted 
curve), respectively. For the three on the right, I = —100, —99 
and -98, respectively. For the two in the middle, I = and 2, 
respectively. 



We have varied the perturbation strength, with cr from 
0.01 to 5, and found qualitatively similar results. These 
results show that for not long times, the two Fock states 
\L) and | — L) are more stable than other Fock states 
with \l\ not close to L. 

There is also some difference between the fidelity decay 
of the two initial states \ — L > and \L >, as shown in 
Fig.Hl that is, the fidelity of \L) decays more slowly than 
that of I — L). This can be understood from the form 
of the Hamiltonian in Eq. (|4|). Indeed, for the state \L), 
the first two terms on the right-hand side of Eq. ^ give 
{L + gcL), while for | — L) they give {gcL — L). Since 
L + gcL > \gcL — L\, the state \L) is less perturbed than 
the state | — L) for the same change of the parameter K. 

We have also studied fidelity decay of Fock states 
|/) with \l\ close to L. Figure [2] gives the cases of 
I = 100,99,98, and 97, which shows that the rate of 
the initial decay of fidelity increases with decreasing /, 
and the oscillation of M{t) appears even at / = 99. For I 
close to —L, the situation is similar. 

In Fig. [31 we show fidelity decay at different values of 
K corresponding to the same fixed time t. For a = 0.01 
shown in the upper panel, at t = 1000 the fidelity of 
the two Fock states \L) and | — L) is much higher than 
that of the other three Fock states for K between 0.5 
and 3. With increasing the perturbation strength cr, the 
fidelity at this time becomes smaller. For a = 0.04 shown 
in the lower panel, M(IOOO) of | — L) is already small, 
while Af (1000) of \L) is still high within some windows 
of K. This is because the fidelity of | — L) decays faster 
than that of \L) (see Fig. [1] and discussions in a previous 
paragraph) . 

In quantum regular systems, the fidelity of initial co- 



where z = 



'cot(6'/2). This gives 



|2(/+L) 



(l + |zP) 



(2L)! 



{L + iy.{L-iy. 



(9) 



Considering the limit \z\ — )■ oo, it is easy to check that 
\L) is also a coherent state. Other Fock states are not 
coherent states; some examples are shown in Fig. |4l 

One may expect that, since the two Fock states \L) and 
I —L) are coherent states, the initial slow Gaussian decay 
of their fidelity, as shown in Fig. [1] might be explained 
quantitatively by making use to results given in Ref . . 
However, detailed analysis show that the situation here 
is more complex. Indeed, in the derivation of Gaussian 
decay given in Ref. 18], a property of coherent states is 
made use of, namely, a Gaussian form of the expansion of 
a coherent state in the basis of Fock states. This property 
is possessed by most of the coherent states, but the two 
states \L) and | — L) are themselves Fock states. If one 
tries to use a 5-function for the expansion, going along 
the line of Ref. [l^ , it is found that their fidelity has a 
decay rate which is equal to zero. Hence, this approach 
can not explain quantitatively the initial Gaussian decay 
of the fidelity of the two states \L) and | — L). But, it 
indeed gives a qualitative explanation to the fact that the 
decay is slow. 



IV. FIDELITY FOR DIFFERENT INITIAL 
STATES 

In practical situations, experimentally prepared initial 
states are usually not exactly the same as the expected 
ones. Hence, in addition to perturbation, one should also 
consider small change in the initial state in the study of 
fidelity, i.e., considering fidelity of a more general form 
than that in given Eq. ([T]). In this section, we consider 
such a more general fidelity and use it in the study of the 
stability of Fock states. 
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FIG. 5: Variation of \mik{t)f with t for « = -31. It has 
relatively regular peaks when \k\ is close to L. Parameters are 
L = 100, K = 2, Qc = 0.17, for which the classical counterpart 
has regular motion, and a = 0.5. 



A. Fidelity for different initial Fock states 

Suppose the expected initial state is a Fock state \k), 
while what is really prepared is the state 

l*o) = EcHO (10) 

with \Ck\ close to but smaller than one. The expected 
state at time t is the time evolution of |fc) under Hq, 
while the real state is the time evolution of \'^o) under 
H. Hence, one should consider the following; more general 
form of the fidelity amplitude (see Ref. [13|), 

fit) = (*o|e'^*e-^^"*|fc) =Y.^tmk{t), 

I 

where 

mikit) = (/|e'^'e-'^°*|fc). (11) 

The quantity mik{t) can be regarded as a generalized 
echo, from an initial state |fc) to a final state \l). Its 
absolute value square, Mik{t) = |m;fc(t)P gives the prob- 
ability for the final state to be found in |^), if the initial 
state is \k) and the dynamics is governed by Hq for the 
first time interval t and by (—H) for the second interval 
with the same length. For k = I, mik{t) is just the fi- 
delity amplitude in Eq. ([1]). In this section, we are more 
interested in the case of fc 7^ Z, for which |m;fc(0)p = 0. 

For a quantum regular system Hq, Mik{t) oi k ^ I 
may be considerably large within some time windows, 
due to the peculiarity of integrability. [If Hq is a quantum 
chaotic system, Mik(t) usually approaches its saturation 
value soon and fluctuates around the saturation value.] 

We would like to mention that the quantity mik (t) also 
appears, when one considers the difference between the 



FIG. 6: Same as in FigElfor I = -L = -100. 

expectation values of an arbitrary observable A for the 
same initial state under two slightly different Hamiltoni- 
ans. Indeed, consider an initial state |fc) and the expec- 
tation values of A in the two systems 

A^°it) = {k\e^"«'Ae''"'''\k), (12) 
Af,{t)^{k\e^"'Ae~^"'\k). (13) 

Inserting the identity operator e~'^°*|Z)(^|e*^''* in 
Eq. (fT3|) before and after the operator A, we obtain 

w 

where the prime over the sum implies that / = k and 
V = k can not hold at the same time and 

Afi?{t) = {l\e'"'>' Ae-'"'>'\1') (15) 

is a quantity given in the system H{). 

B. Numerical investigation 

Now we study the behavior of the quantity M/^ [t) in 
the two-component BEG system discussed above. As 
shown in Sec. IIIIl for some initial times the two Fock 
states \L) and | — L) are more stable than other Fock 
states 1^) with |Z| not close to L, in the sense that the for- 
mer's fidelity is higher than the latter's. In this section, 
we show that our numerical investigations in the quantity 
Mik{t) give consistent results. Specifically, Mik{t) with 
either \l\ or |fc| close to L behave more regularly than 
those with both |^| and |fc| far from L. 

Let us first consider Mik{t) with |/| far from L. An ex- 
ample is given in Fig. [51 which shows variation of Mik{t) 
with time t for L — 100, / = —31, and k from -100 to 
99. It shows that for k close to —L, Mik{t) has quite 
regular peaks, and with increasing k from —L to 0, the 
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FIG. 7: Variation with k for the times corresponding to the 
centers of the first and second peaks of \mik{t)\'^ shown in 
Fig.El 



FIG. 8: Values of Sk = 'Lk'<k Mik'{t) for I = -L. Param- 
eters are the same as in FigO It shows that at each of the 
times, Mik{t) is concentrated in a small region of k. 



peaks become more and more irregular. Similarly, for k 
decreasing from L to 0, the peaks of Mik{t) also becomes 
more and more irregular. 

On the other hand, for I — — L — —100, Mik{t) behaves 
regularly for all the values of fc, as shown in Fig. [51 Here, 
the basic feature is that for each value of fc, Mik{t) is con- 
siderably large only within some time intervals. Specif- 
ically, for k = —100, Miic{t) = 1 at t = as a trivial 
result of fc = I, then, it decays and remains small until 
t ~ 1370 after which there is a revival, forming a second 
peak centered a.t t ~ 1450. The second peak of Mik{t) 
oi k = —100 splits into two peaks when k = —98, hence, 
Mikit) has three peaks for k = —98. Interestingly, with 
further increasing fc, the first peak of Mik{t) moves to the 
right while the second peak moves to the left. The two 
peaks meet at a value of fc a little smaller than 79. In 
Fig. [3 we show variation of the times corresponding to 
the centers of the first and second peaks with increasing 
k. 

The structure of the peaks shown in Fig. [S] suggests 
that at each fixed time t, \Mik{t)\ of I — —L may be 
concentrated in a relatively small region of k. Indeed, 
this is confirmed by our further numerical simulations. 
We have calculated the quantity 

= ^ Afifc,(t), (16) 

k'<k 

which gives the total probability for k' < k. Examples 
for some fixed times are given in Fig. [51 which shows that 
|M;fe(i)| of I = —L is indeed concentrated in a relatively 
narrow region of k for each time. 



V. CONCLUSIONS 

By numerical simulations, we have studied quantum 
Loschmidt echo or fidelity decay of initial Fock states in 
a two-component BEG system, whose classical counter- 
part has regular motion. Our results show that, for some 
initial times, initial Fock states with all the atoms in one 
component of the BEG are more stable than Fock states 
with atoms distributed in the two components. This im- 
plies that one-component BEG might be more stable than 
two-component BEG. 

We have further investigated this issue by considering a 
more general form of the fidelity, i.e., fidelity for different 
initial states. Numerical computations of the general fi- 
delity show consistent results, namely, initial Fock states 
with all the atoms in the same component behave more 
regularly than other Fock states. 
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APPENDIX A: AN EXPERIMENTAL SCHEME 
FOR MEASURING FIDELITY 

In this appendix, we discuss in detail an experimental 
scheme for measuring fidelity decay in a BEC system, 
which is briefly sketched in the last section of Ref. Q, 
and give an explicit expression of the fidelity in terms of 
measurable quantities. 

Experimental schemes for measuring fidelity have been 
discussed by several groups and basically three types of 
schemes have been proposed [1^. In the first type of 
scheme, a quantum system is considered, which is com- 
posed of two subsystems (or has two degrees of freedom) 
pgl [30I [39I lioj . The system is assumed to have a time 
evolution such that the fidelity of the first subsystem is 
given by the reduced density matrix of the second subsys- 
tem [41I, |43]. Then, measuring properties of the second 
subsystem which can be small, fidelity of the first subsys- 
tem which may be large can be obtained. This scheme is 
adopted in the experiments in Ref. [H, . To be spe- 
cific, one may consider a Hilbert space which is the direct 
product of a two-dimensional subspace with basis states 
|1) and 1 2) and a second subspace for state vectors 
The Hamihonian has the form H = + i?2|2)(2|, 

where Hi and H2 act in the second subspace only. For 
an initial state |*(0)) = (|1) + \2))\ipo)/V^, the state 



at time t is \^{t)) 



1 

V2 



-iHit 



l^o)|l> 



-iH2t 



l^o)|2)) 



Then, the fidehty amphtude ('0o| 



1^0) 



the 



second subspace is related to an off-diagonal element of 
the reduced density matrix in the first subspace. 

In the second type of scheme, a special kind of system 
is considered, for which the time evolution can be such 
controlled that the system evolves under Hq for a time 
period t, then, evolves under {—Hq — eV) for a second 
period t. The fidelity is then just the survival probabil- 
ity, i.e., the probability for the final state to be found in 
the initial state. In the third type, classical waves are 
employed, which evolves according to a dynamical law 
mathematically equivalent to Schrodinger equation [45| . 

Now we discuss the scheme briefly mentioned in 
Ref. We propose to use a setup similar to that used in 
Ref. Hg. Consider a BEC (e.g. which is optically 

cooled and trapped, then, transferred into a double-well 
potential. The double-well potential can be created by 
deforming a single- well optical trap into a double- well po- 
tential with linearly increasing the frequency difference 
between the rf signals [46| . Near- resonant coupling fields 
are applied to the BEC in the two wells with slight dif- 
ference in strength. Finally, simultaneously switching off 
all the external fields and letting the BEC expand freely, 
interference pattern of the BEC can be observed. The 
wells should be deep, such that the total density remains 
approximately a constant, and the atom numbers in the 
two wells are required nearly equal. 

At the initial time to, suppose the state of the sys- 
tem is a product state \^{to)) = \(j){to))\^fj{to))\^p(to)) , 
where \<p{tQ)) is the internal state of the atoms, e.g., with 
all the atoms in the same hyperfine internal state, |V'(to)) 
describes the motion of the center-of-mass degrees of free- 
dom of the atoms, and \^p{to)) represents the field form- 
ing the optical trap. We assume that the field of the 



optical trap is not entangled with the BEC in the exper- 
imental process and shall omit the term \^p{t)). 

From time to to ti, the potential of the optical trap 
is deformed into a double-well potential. If the internal 
state of the atoms is not influenced in this process, at t = 
ti, one has = |0(ti)) [|Vi(Ri, ti)) + \M^-2, h))] , 

where Ri and R2 indicate spatial locations of the two 
wells, respectively. 

From time ti to t2, near- resonant coupling fields can 
be applied to the condensates to couple the two hyperfine 
states. The coupling fields have a slight difference in 
strength in the two wells. We assume that the near- 
resonant coupling fields can be treated as classical fields 
and do not induce tunnelling between the two wells. The 
internal states of the condensate in the two wells will then 
evolve differently. Thus, for t E (ti,t2), 

\^{t)) = \MmMRl,t)) + |</>2(t))|V'2(R2,t)). (Al) 

In the case that the internal degrees of freedom are not 
coupled to the center-of-mass degrees of freedom, \(j)j{t)) 
has unitary time evolution 



|0,(t)) =[/,(ti,t)|<^(ti)). 



(A2) 



with J = 1, 2 indicating the two wells. The internal state 
can be expanded in the Fock states \l), 



(A3) 



At t = t2, one can simultaneously switch off all the 
external fields and let the two clouds of BEC expand 
freely. For t > t2, Eqs. IXTjl and ((X3|) are still valid. 
Substituting Eq. (|X3)) into Eq. (KA\ . one has 

= E [d'i'\t)\MR,,t)) + 4'\t)\MR2,t))] |0.(A4) 

I 

Suppose the single particle states for |?/'i(Ri,t)) and 
|V'2(R2,t)) are xi(x, t) and X2(x, t), respectively. Then, 
the probability of finding a particle at a position x is 



P(x, t) = IxiP + IX2I' + 2Re [/(t)xiX2 
where m^j:,d['\t)d'i'^*{t). 



(A5) 



Making use of Eqs. (|A2[) and (|A3p . it is seen that 

fit) - (0(ti)|{/t(ti,t)t/i(ti,t)|0(ti)), (A6) 

which is a fidelity amplitude. Since there is no coupling 
field beyond t2, Ui{t2,t) — U2{t2,t) and, as a result, 
/(t) = /(t2) for t > t2. Then, Eq. (jM]) can be written 
as 



2|/(t2)|Re 



=i«/(t2) 



X1X2 



P(x,t)-|xip-|X2p, (A7) 



where ^^/(^2) is the phase of /(t2)- Therefore, the value 
of |/(t2)| can be obtained by measuring the interference 
pattern of the two expanding clouds of BEC. 
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